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The Wigner function is known to evolve classically under the exclusive action of a quadratic 
Hamiltonian. If the system also interacts with the environment through Lindblad operators that are 
complex linear functions of position and momentum, then the general evolution is the convolution 
of a nonhamiltonian classical propagation of the Wigner function with a phase space Gaussian 
that broadens in time. We analyze the consequences of this in the three generic cases of elliptic, 
hyperbolic and parabolic Hamiltonians. The Wigner function always becomes positive in a definite 
time, which does not depend on the initial pure state. We observe the influence of classical dynamics 
' and dissipation upon this threshold. We also derive an exact formula for the evolving linear entropy 

I as the average of a narrowing Gaussian taken over a probability distribution that depends only on 

. the initial state. This leads to a long time asymptotic formula for the growth of linear entropy. We 

04 ' finally discuss the possibility of recovering the initial state. 

^ ! 
Ph. 

<^ ; I. INTRODUCTION 

i> : 

, The correspondence between classical and quantum mechanics of closed dynamical systems is most perfect for 
quadratic Hamiltonians. In this case, the classical evolution is linear, like its quantum counterpart, and generates an 
' orbit within the group of symplectic (linear, canonical) transformations in phase space [1]. These are directly linked to 
^ the corresponding quantum metaplectic group [2]. Indeed, the evolution operator in any of the usual representations 
is merely the complex exponential of the classical generating function [3]. Of course, quadratic Hamiltonians are a 
very special case, but they include the ubiquitous harmonic oscillator, the parabolic potential barrier and the free 
■ particle, which form adequate starting points for the analysis of more complex motion. 
. The Weyl representation of an arbitrary quantum operator A is 

m ■ 

^(x)^ y"dg'(g+^|l|g-^)exp(^-z^), (1) 

^ . that is, A is represented in phase space, x = (p, g), by the Weyl symbol A(x). The Wigner function VF(x) is then 
^ ' the Weyl symbol for p/2-Kh, where p is the density operator. Like all Weyl symbols, the Wigner function propagates 
^ , classically under the action of a quadratic Hamiltonian [2] [3] : 

-W^t(x) = {ff(x),VFi(x)}, (2) 

introducing the classical Poisson bracket on the right side [1] and i/(x) = x • Hx, the Weyl symbol of the quadratic 
Hamiltonian. The symbol • stands for the inner scalar product. Hence one has Wti'x.) = Wo(R— tx), where 



00 

o 



Rt = exp (2JHt) (3) 

is the 2x2 matrix giving the classical Hamiltonian time evolution of a phase space point x. Actually, this propagation 
is also shared by the Fourier transform of Wt (x) , 

^(0 = A exp f-k A xVw, (4) 



27rn ./ \ h 

i.e. it is also true that 

d 



Q-Wt{i)^[H{0,Wt{i)]. (5) 
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Note that -ff(^) must be read literally as the classical Hamiltonian -ff(x) taken at the point ^ = {£,p,S,q), which is in 
general different from the chord transform H{^) of H{x). In other words, one has also Wt{^) = Wo(R_t$). Above we 
have made use of the wedge product, 

CA^ = Cpq-^gP= (JO^x = Je • X, (6) 

also defining the transpose of a vector (.)^ and the matrix J. The semiclassical background for the symplectic 
invariance of both the Wigner function and its Fourier transform is that the Weyl phase space coordinate x may 
be associated to pairs of points in phase space, x±, by x = (x+ + x_)/2. The conjugate variable to this centre is 
the chord ^ = x_|_ — x_. The linear motion of both the chord ^ and the center x is the same as for each individual 
phase space point x+ or x_. We will refer to W{^) as the chord function as in [4], though it is also known as the 
characteristic function in quantum optics. 

The question that we address is to what extent can the simplicity and generality of symplectic motion of closed 
quantum systems be incorporated within the description of systems whose coupling to the environment cannot be 
ignored. In this case the evolution is no longer unitary, unless the full Hamiltonian of the system combined with the 
environment is taken into account. All the same, a certain measure of generality is restored by the assumption that 
the density operator is governed by a Markovian master equation [5], 

f = -^[^'^] - 2^^ 2L,pLt-L]L,p-pLtL,. (7) 
j 

If the further assumption is made that each Lindblad operator, Lj, is a linear function of p and q, we will show that 

the evolution of Wt(^) is the product of the classically evolved Wo(^) with gt{C), a Gaussian centred on ^ = that 
has diminishing width. One can then generalize in a straightforward way the exact solution given by Agarwal [6] for 
the Wigner function. This is a convolution of the Fourier transform of gt{C) with the classically evolved Wo(x). In 
other words the Wigner function is coarse-grained by a widening Gaussian window. 

A simple example of symplectic evolution of an open quantum system is that of a dust particle interacting with 
air molecules, or radiation, so that in the absence of gravity, H = fP' /2m and the interaction with the environment 
depends basically on the particle's position: L = rjq, where r? is the coupling parameter. This example is discussed by 
Giulini et al. in [7]. Another important example is that of an optical field, say an arbitrary superposition of coherent 
states interacting with thermal photons. In terms of real variables, the internal Hamiltonian is just H = uj{p^ -\- q^)/2, 
i.e. the harmonic oscillator. The Lindblad operators in this case are known to be 7(n + l)a/2 and 7na1'/2, where 
and a are the usual creation and annihilation operators, n is the average number of thermal photons at the frequency 
oj of the cavity mode at temperature T and 7 is the decay rate [8]. 

Recently Diosi and Kiefer [9] (hereafter labeled DK) showed in the case of the first example that the Wigner function 
of any pure state becomes positive within a definite time. Thus the Markovian interaction with the environment has 
the effect of erasing the interference fringes characteristic of quantum coherence and from then on the effect of coarse- 
graining on the Wigner function is not distinguishable from that of a classical Liouville distribution. How general 
is the DK scenario? In section // we present the exact solution of the Lindblad equation for general quadratic H 
and arbitrary complex linear Lindblad operators Lj = Xjq + fJ-jp. In section /// we explicit the underlying classical 
structure of the solution. Then, in section IV, we use the properties of the convolution to generalize the positivity 
time of DK in the case of arbitrary quadratic Hamiltonian and nonhermitian Lindblad operators. Furthermore we 
make the much stronger statement that the Wigner function cannot be positive before this threshold, unless the 
initial distribution is a Gaussian. We discuss the consequences of this statement through the example of a bath of 
photons. In section V we specify the behaviour of the positivity threshold for each type of quadratic Hamiltonian. 
It turns out that a nonzero dissipation coefficient implies that positivity is reached exponentially fast. The positivity 
threshold is in general inversely proportional to the dissipativc coefficient. However it becomes inversely proportional 
to the Lyapunov exponent if the latter is greater than this coefficient, that is in the case of a hyperbolic system 
(i.e. the inverted oscillator) in the weak coupling limit. Then positivity can be reached much faster than in the 
corresponding elliptic case (i.e. the harmonic oscillator with the same coupling constants). Though all the formulae 
presented here are appropriate for a single degree of freedom, the generalization to higher dimensions is discussed in 
this section. In section VI we derive a general formula for the growth of the linear entropy (1 — Tr pt^), with respect 
to the initial density operators, for each choice of the quadratic master equation. We also show that for long times 
the growth of linear entropy attains a universal form. Finally in section VII we point out that this general solution is 
obviously reversible, giving a very synthetic inversion formula which generalizes previous work about quantum state 
reconstruction [10]. 
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The generalization of the convolution as exact solution of the master equation (7) when H is not quadratic, or for 
nonlinear Lj, is not obvious. However, the approximate semiclassical theory developed by one of the present authors 
[11] has no such constraint. Its compatibility with the present theory is the subject of a companion paper [12]. A 
simpler version of the present work, for the restricted case of hermitian Lj, can be accessed in [13]. 



II. EXACT SOLUTION IN THE QUADRATIC CASE WITH DISSIPATION 



We derive here the exact solution of the Lindblad equation in the case where the Hamiltonian is quadratic and the 

Lindblad operators arc complex linear forms in q and p. 

Taking the Weyl-Wigner transform of equation (7), i.e. associating the Weyl symbol A{x.) to each operator A, as 
defined in (1), and using the product rules [14] for operators, we obtain. 



dWt 
dt 



(x) = {j?(x),m(x)} + ^(ji;M;.) 



dWt 

x.^(x)+2W^,(x) 



+ 



d'^Wt 



(x) Jl' + Jl" 



ax2 



(x) Jl" 



(8) 



Here the Lj (x) = 1^ • x + il" ■ x are the Weyl symbols of the nonhermitian linear Lindblad operators. We use the 
notation 



X = 1^1 and 1' 



A' 



1" 



A" 



(9) 



It is well known that in this case H{x) and -Lj(x) can be identified as the classical variables corresponding to H and 

Lj. Note that if I" = then the second term in (8) cancels. It will become clear that this term is responsible for 
dissipation in the evolution of W((x) and we define the dissipation coefficient, a = X^j(Jlj • lj), which is zero in [13]. 

It is actually easier to solve the evolution equation for the chord function H^(^), 



dWt 

dt 



(0 = HaM^.(0}-«e-^(?)-|:E 



(10) 



as derived in the appendix A. We guess a solution of the form 



dt' 



(11) 



where is a linear evolution of ^, which will be explicited a posteriori, such that 

^o = e 



(12) 



Then, inserting the form (11) of Wt in equation (10) and dividing both sides by the exponential of (11) leads to the 
following left part. 



dWo 



- E / [2 (i;- • ^t'-t) I'j ■ {-it'-t) + 2 (ij • 6'-t) 1;' • (it'-t) 



dt', 



(13) 



which must equal the right part 



^ - '^^^-o-4Efe-^)^ + (i"-^)1^^o(^-o-4^o(^-oE 



f 2 [(1;. • 6'-t) 1 • • (-2JH6'-t - a^t'-t) + (1; • it'-t) i;- • (-2JH6'-t - a^t'-t) 
Jo I 



2h 

dt'. 



(14) 
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We have used 



and 



(wo(6)) = 6-^(6), (15) 



f ~ -I dWt 

p(0,^^t(o} = -2JH^^(e). (16) 

Hence the ansatz (11) is a solution of (10) if fulfills 

6 = 2JHCt + a6- (17) 

Thus, we can write explicitly, 

6 = e"*Rtf , (18) 

where Rj, defined in (3), gives the purely Hamiltonian evolution and the dissipation term a leads to a classical 
nonhamiltonian expansion (a > 0) or contraction (a < 0) of the chord variable ^. One should be aware that although 
the Hamiltonian part of the evolution of ^ is shared with the one of the phase space point x, the effect of dissipation 
is inverted, as it will be explicited soon. 

The argument of the exponential in (11) is a quadratic form in ^, so the solution can be written 

mo = M^o(C-t)exp (^-^i ■ , (19) 

with M(t) a real, time dependent 2x2 matrix, which can naturally be decomposed into 

M{t) = ^ M,(t) = ^ dt'e2"(*'-*)R?;_,l,ljR,,_„ (20) 

so that each Lindblad operator contributes a Gaussian to (19). 

Now, back into the Weyl-Wigner representation by using (4), one obtains the solution of (8), 

where we have defined 

Mj{t) = -3M{t)J 
Mj(t)-i = -JM(i)-^J, (22) 

with the symplectic matrix J defined in (6). Wo{x) is the initial Wigner function and the convolution Gaussian turns 
into a Dirac 6 function as t goes to 0. We have equivalently, 

Wt(x) =t«t(x_t) (23) 

1 r / I 



: I Wo (y) exp (-|:(y - x_t) • M{t)-\y - x.*)) dy, (24) 



2nhy/detMj{t) 
with 

xt = e-«*RtX (25) 

and 

M{t) = -e2"*R!;tJM(t)JR_t = -Mj{-t). (26) 

Hence the solution is a convolution with a Gaussian which broadens in time, composed with a backwards nonhamil- 
tonian evolution of the phase space variable x. As mentioned earlier, the Hamiltonian part of this classical evolution 
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is the same as in the chord space, whereas the dissipative part has the opposite effect: dissipation (a > 0) will shrink 
the phase space variable, and thus expand the distribution Wt- 

This solution, which has been derived in the case of an homogeneous quadratic Hamiltonian, can be generalized 
easily to a quadratic Hamiltonian with a linear part. One then has to be aware that the matrix Rj, which appears in 
the exponential Lindbladian damping, strictly corresponds to the classical motion of the chord, which is determined 
by the homogeneous part of the Hamiltonian H. This remark is important, for instance, in the parabolic case, say a 
particle with a linear potential, where the motion of the chord disregards the potential. 

Since | detM(f)| grows with time, one can conclude, following DK [9], that the solution (21) becomes positive after 
a certain time. Indeed it is, after rescaling the variable, the convolution of the initial Wigner function with a Gaussian 
of broadening size, which smoothes out oscillations around zero. It is the property of symplectic invariance of the 
Wigner function that DK employ to prove strict positivity in a specific case that is now extended to its broader 
context. Moreover we shall give in section IV the much stronger resiilt that the Wigner function cannot be positive 
before the DK time, which does not depend on the initial pure state, unless it is a coherent state. We then give in 
section V the general behaviour of this positivity threshold for different dynamics, namely when H is elliptic, as for 
the harmonic oscillator, hyperbolic, as for the scattered particle, or in the parabolic intermediate case which includes 
the system studied in DK. The next section explicits the formal correspondence between this problem and a classical 
Brownian motion described by a Langevin equation, as also did Agarwal [6] for his solution. 



III. CLASSICAL CORRESPONDENCE 



Since Eq. (8) is a Fokker-Planck equation, it can be interpreted as the evolution equation for the probability 
distribution of a classical Brownian motion defined by a Langevin equation. This correspondence gives a simple 
classical interpretation of the problem. Hence the decoherence may be seen as a diffusion induced by some random 
force, and dissipation can be interpreted as a classical viscosity, although it is always accompanied by another diffusive 
term. The only feature which cannot be assigned a classical meaning is the Wigner function itself, which, as a pseudo- 
probability distribution, can have negative values. 

One can check easily, see for instance [15], that the following Langevin equation. 



dq 



(27) 



induces equation (8) as a Fokker-Planck counterpart. The "Brownian forces" fm{t) and gm{t) verify 



{fmit')fnit")) = d^.nS{t) 

{9m{t')gn{t")) = Sm,nS{t) . (28) 
{fm{t')gn{t")) = 

The f,n{t) correspond to the diffusion induced by the nondissipative real part of the Lindblad operators, whereas the 
gm(t) correspond to the diffusion induced by dissipation. It can easily be verific;d that the Fokker-Planck equation is 
symplectically invariant, so one is allowed to perform the following change of coordinates. 



p = p-mr.\ (29) 

q = q 

to the above Langevin equation, which then turns into the following more intuitive form, where the dissipation depends 
only on the momentum: 

Here, the matrix for the transformed Hamiltonian H is 

Hii Hi2 

TT TT I a^+4QHi 
Hi2 H22 + 4Hii 

^'m — X'm~ A*m ^^'^ respectively for A^, and a = 2a. 



(31) 
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From this classical picture we can interpret the general behaviour of the solution of the Lindblad equation. In the 
case of a closed system, remember that the Wigncr function undergoes a Liouville, unitary, propagation. Now the 
system is coupled to an environment, i.e. there arc nonzero Lindblad operators. If dissipation of energy is neglected, 
these operators are Hermitian, so there is no imaginary part of the vectors 1. Then the effect of the environment 
over the system can be interpreted as a diffusion process corresponding to random forces in the Langevin equation. 
Formally, it corresponds to the initial Wigner function being convoluted with a Gaussian which broadens with time. 
If one now takes into account the dissipation induced by the environment, allowing the Lindblad operator to be non- 
hermitian, a viscous term appears in (28) and (30), meaning that the classical trajectories on which the distribution 
travels arc drifted to lesser or higher energy, according to the sign, + or — , of the "viscosity" a. Indeed, the dissipative 
linear motion governed by the nonrandom terms of (28) is just that of (25). One should be aware that it is only a 
formal classical scheme, and that the viscous term might have a purely quantum origin. For instance, in the case 
of photons in a cavity with dissipation, whose Lindblad operators are explicited in section V — A, the viscosity is 
related to spontaneous emission, which breaks the symmetry between emission and absorption. Then the classical 
trajectories of the above equation spiral in towards the origin, although a semiclassical theory would lead to no 
dissipation. The opposite case would be an amplified cavity, where the trajectories would spiral out. Notice that 
this viscous term always goes along with a supplement of diffusion, which can be interpreted as a consequence of the 
fluctuation-dissipation "theorem" [16]. 



IV. UNIQUE POSITIVITY TIME FOR ANY INITIAL STATE 

We have seen in section II that, because of the Lindbladian part of the master equation, the pure state Wigner 
function is convoluted with a Gaussian whose width grows with time. It has been pointed out by DK [9] that at 
the time tp at which the width of the Gaussian reaches h, that is when it becomes the Weyl representation of some 
coherent or squeezed state, then its convolution with the initial Wigner function, Wo(x), is a Husimi function [17] [2] 
[4] of the initial state, or a Q function in the language of quantum optics. It is a well known property that the Husimi 
function is positive, so we have Wt^ > 0, and obviously, since the Gaussian is strictly broadening, Wt > for t > tp. 
It has already been emphasized by Leonhardt et al. [18] that the form of the Wigner function after interaction with 
a dissipative environment can be read as an intermediate phase space distribution W{x,t = 0, s), with an s which 
depends on the dissipation rate. Thus s — corresponds to the (initial) Wigner function and s = — 1 to the Husimi 
function of DK (in [18] the role of the environment is played by the imperfections of a beam separator) . 

We shall now prove that the Wigner function can never be positive before the positivity threshold tp, unless it is 
positive from the beginning, that is, unless the initial state is a Gaussian state. Indeed, if an initial pure state \ipo) is 
not a Gaussian, then it was shown by Tatarskii [19] that the initial Wigncr function 14^0 (x) has negative parts. But 
it is also true that nongaussian Husimi functions necessarily have zeroes, as shown in the appendix B. Since Wt^ is, 
up to a symplectic transform, a Husimi function, there exists xq such that 

W^t,(xo)-0. (32) 

Let us now investigate the case t < tp. Then Wt{x) given by (24) is a convolution of Wo{x_t) with a Gaussian of 

width smaller than h. The point is that one can then convolute again with exp ^— ^x_( • M(t)~^x_t^ , with the real 

parameter k chosen so that the output is also a Husimi function (3(x_t). To show this we refer to the simple general 
relation: 



(33) 



/ dy exp(-(x - y) • M(x - y)) exp (-y • AtMy) = ^ exp ( -x • — Mx ). 

J (H-K)VdetM V l + « / 

This Husimi function can be identified with Wt^ through a symplectic transform x_t — > x', so we have 

g(x^) = Jdy wti^'o - y)exp (-^y • M(t)-V) = 0. (34) 

Now it is obvious from (34) that Wf, hence Wt, must have a negative part for all t <tp. 

Let us emphasize the striking consequence of this result: the positivity time does not depend on the initial distri- 
bution, as long as it is not a Gaussian one. The following example shall illustrate this remarkable property in a more 
transparent way. 
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A. Example 



We start from the familiar superposition of two coherent states, i.e. ground states of the harmonic oscillator, 

displaced to the phase space points = (0,±C): Itpa) = {\0 + \ — C))/ y/{2), in the context of photons in a cavity 
with dissipation. The Wigner function, which is a particular case of Eq. (21), is the sum of three terms; two of these 
correspond to the coherent states taken independently, and the third term comes from their interference: 



with 



and 



Wt(x) = Wc(x)+M^_^(x) + iy,(x), (35) 



H^C(X) = ^ exp (-1,^) exp (-!(,- e-V^C)^ ) (36) 



W.i.) = i|exp (-|(/ + .^)) exp (-2(1 - ^)C^) cos (l^Cp), (37) 

where 7, defined in the introduction, corresponds to 2a, and (3t = 2n(l — exp (— 7t)) + 1. The function is normalised 
hy N = (1 + exp {—Q'^/Kj) . Obviously the minimum values are concentrated on the line q = Q, where the expression 
simplifies into 



Wt{v, g = 0) = ^ exp {-^p" 



exp ( -2(1 - cos — + exp 



fit J V A ' ^ V A 



(38) 



The Wigner function becomes positive when 



that is at tp = l/7log(l + l/(2fi+ 1)), which indeed does not depend on the initial spacing ( of the two coherent 
states. On the other hand, the position (^^,0) of the closest zero, given by the first minimum of the cosine at that 
time, 

gets further away as C becomes smaller. This shows that the negative regions that remain until tp may be so shallow 
as to be practically irrelevant. 

The chord representation reveals how the positivity threshold is related to a more reasonable estimate of the 
decoherence time. The expression of the chord function in our example is 

WoiO = A -P (-Q 2 cos (^-^) + -L exp (-^^^) + A -P f-^^) • (41) 

2nh \ 4hJ \ h J 2TTh \ m J 2-Kh \ m J ^ ' 

The first two terms, distributed around the origin, correspond to the two coherent states taken separately, whereas 
the last two terms, distributed around x.^ and — x^, describe the quantum interference between them. The positivity 
time corresponds to the moment when the Gaussian in (19) damps everything outside a region of size h in the chord 
space. However, both interference terms in this example will be damped much sooner if C, is large enough, indicating 
an overall loss of coherence. 

It should be remarked that one just has to study a specific example, here the twin coherent states, to get the 
positivity time for any nongaussian initial pure state. 



V. BEHAVIOUR OF THE POSITIVITY THRESHOLD 

Positivity is attained when the determinant of the matrix M(<), that is the determinant of M(— t), is equal to 1/4. 
Then the expression of the solution (24) is indeed a Husimi function. The matrix M is defined by (20), so, to calculate 
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it, wc notice that JH, in the expression (3) of Rt, can be diagonahzed in most cases, that is when its two eigenvalues 
are finite and different. We then define the matrix P such that 



2JH = P-^DP, (42) 

with 

and cr = 2\/— det H. Notice that since H is symmetric, then JH has a null trace, and so has D. The dissipation 
parameter, q, and a are basic elements for the description of the evolution of Markovian quadratic open systems. 
Then, using (20) one can easily derive the expression of M, 

Mit) = P^ ij^ " P' (44) 

where A is defined by 

A=(p-i)^5](l^(l^)^ + lJ(i;)^)p-\ (45) 

j 

The corresponding quadratic form actually depends on the classical motion of the chord, determined by the homo- 
geneous part X • Hx of the Hamiltonian. One then has to separate different cases: the elliptic case, det H > 0, the 
parabolic case, detH = 0, and the hyperbolic case, detH < 0. 

Since the Wigner function is symplectically invariant, one just has to treat the simplest expression in each case, 
respectively the harmonic oscillator, (x) = p^/2 + the particle in a hnear potential, -ff (x) = p^/2 + q, and the 
scattered particle, H{x) =pq, which is symplectically equivalent to p^/2 — g^/2 ^. 



A. Harmonic oscillator 



In this case, the Hamiltonian reads 



and the matrix P is 



Then the determinant of M.{—t) reads 



H{k)=u;(^ + ^] (46) 



, N e^«*-2e2«*cos2wi+l , . e^"* - 2e2«* + 1 . . 

det M{-t) = o , 0, An A22 A12 A21 , 48 

4(a^ + u-') 4a^ 

where in this case the matrix A is complex: 
using (9). 



^The reduction of a linear Lindblad operator under a syniplectic transformation, x' = Cx, is especially simple in the Wigner 
or in the chord representation. If, instead of (9), we set I/j(x) = Ij A x, then the invariance is obtained by taking I'j = Clj. Of 
course, one must also use the invariance of the classical Hamiltonian -ff(x') = i?(x). 
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In the dissipativc case, that is for a > 0, this determinant diverges exponentially fast, and positivity is attained 
in a time of the order of 1/|q;|. Let us take for instance the bath of photons, then the coefficients of the Lindblad 
operators are 




and the friction a = 7/2. Then 

detM(-t) = ^— ^^(2n+l)2, (51) 

which equals 1/4 at t = l/7log(l + l/(2n+ 1)), as was previously mentioned. Notice that, although the coarse- 
graining grows for ever, the size of the Wigner distribution reaches a finite limit, for the rescaled function (24) is the 
expression of W{yi) and not VF(x_t). 

If, on the other hand, a < then the determinant reaches its limit in a time also of the order of 1/|q;|. We conjecture 
that this limit has a lower bound greater than 1/4. Here rescaling x_t — > x now implies that Wt(x) spreads with no 
bound. 



B. Scattered particle 

The simplest form of the Hamiltonian for a particle scattered by a parabolic barrier is 

fl'(x) = ojpq. (52) 

Then the matrix P is just identity and the matrix A is real, so the determinant reads 

e^"' - 'h '-"' clrl^'l + 1 e''"' - '^('-"' + 1 

det M{-t) = "„ " AnA22 ^ A12A21, (53) 

4(a^ — w^j 4a^ 

As long as a > —u) it grows exponentially, so positivity is always reached. When a < —uj the determinant has again 
a finite asymptotic value. Then the positivity threshold is of the order of l/(|a| — w), which is greater than the 
corresponding elliptic case, with identical Lindblad operators. 

The main difference with the elliptic case appears in the weak coupling limit \a\ <^ u. Whereas positivity of the 
elliptic system will then be attained in a time which is still inversely proportional to the coupling with the environment, 
the positivity threshold of the hyperbolic system will saturate at We conjecture that this will also be the case 
in a more general chaotic system. 

C. Particle in a linear potential 

We now study the intermediate case, which can be represented by the Hamiltonian 

il(x) = ^+g. (54) 

This degenerate case docs not follow our general form for the matrix M, so one has to treat it separately, taking care 
of the linear term (cf. remark of section //). One should remember here that the motion of the chord is given by the 
free motion of the particle, which corresponds to 

R.=(;;). (55) 

Then the damping matrix is 

mt) = ty^ ,^9c.t /f^-.r^ i, (56) 
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where the qI''-' arc polynomials of degree dint and of degree 2 in the coupling constants, say (A', A", /i', p,''). Let us 



take for instance one Lindblad operator with (1')^ = (0 



det M(-t) = - 
4 



AeDt 



1 



1 



4(£)/')2 



„2eDt 



D') and (1")^ 
D' 



D" 



= (e/D",0), with e = 

1 ,A . 1 



±1. Then one has 



2(Li")2 



+ 2 +1 



(57) 



with D = \J D'D" . The limit is always greater than 1/4, with the usual exponential contrast between the dissipative 
case (e — 1) and the excited case (e = —1). In [9], DK study this example with no dissipation, and they find a 
positivity time of the order of I/a/D'. On the following table, one can read different values of the positivity threshold 
for D' = 2, and check that the limit D" ^ is attained gradually: 



D" 





0.1 


1 


10 


100 


e= -1 


0.930 


0.640 


0.244 


0.077 


0.022 


e= 1 


0.930 


1.040 


1.025 


0.752 


0.400 



(58) 



However for large values of D", the positivity threshold will behave like 1/ ^/ D' D" . 

So far the discussion has been restricted to the case of a single degree of freedom. Besides trivial changes of 
factors of 27r/i, the basic form of the solutions of the Lindblad equation in the chord representation (19) and for 
the Wigner function (21) remain unchanged in the case of n degrees of freedom. The evolution matrix now has 
the dimension (2n) x (2n), but it can again be simplified by symplectic transformations. It will often decompose 
into blocks corresponding to the above examples. If every Lindblad vector 1^ is defined for a single block, then its 
contribution to M(t) will be of the same form as that for a single degree of freedom, but otherwise each case must 
be examined separately. Also, 4— dimensional blocks may arise, corresponding to hyperbolic spiral motion, as well as 
singular cases analyzed by Arnold [1]. 



VI. THE GROWTH OF LINEAR ENTROPY 

Besides considering the positivity of the Wigner function, we can use the exact solution (19) to investigate the 
growth of linear entropy 



St = l-Trpt'. 

This is only zero for a pure state, just as for the Von Neumann entropy 

§t = -Trpflnpt. 



(59) 



(60) 



Since the solution is simpler in the chord representation, we make use of the following relations^ for operators A, 
B,. . . represented by chord functions A{^), B{^) . . . 



A ^ A{0, 



(61) 



(62) 



TrA = A{^ = 0), 

TrAB = ^JdiAm-i). 

Therefore, using P =p and pt{C) = 27r?i Wt(^), we obtain 

Trpt = 2TThWt{^ = Q) = l, 



(63) 
(64) 

(65) 



^Note that there is a misprint in formula (6.24) of [4] 
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and 



Trpt- 



= 2-Kn j 



(66) 



In the case of the solution (19) of the Lindblad equation, we thus have 



Woi^ = 0) = — 



and 



1 



(67) 



(68) 



if the initial state is pure. 

After changing the variables of integration in (66), replacing in (19) by the initial chord we obtain 



W^o(OI'exp(ir-M(-t)r 



(69) 



where M(— i), defined in (20), is negative definite. 

Thus Tr pt^ is just the average of a rcscalcd Gaussian whose width narrows in time in a manner that depends 
exclusively on the particular form of H{x.) and the linear Lindblad operators Lj{'x.). The initial state merely determines 
the probability density employed in the calculation of the average. A general asymptotic behaviour can be predicted 
for this formula, because the width of the Gaussian generally shrinks as t grows. If the contraction is sufficient, the 
expression (69) will tend to 



Tr p^ ~ e^"* 



d(, exp 



irhe 



2at 



VdetM(-t)' 



(70) 



by using (67). It can be explicited, using (43) and (44), in terms of the dissipation, a, and the basic Hamiltonian 

parameter, a: 



Trpt' 



-4Qt _n^-2at I 



p2crt_i_p — 2tTt ' 



(71) 



+1 



4(a^+<7^) 



All A. 



-4at_o^-2af 



11^22 



AioA. 



12^21 



If the underlying classical system is elliptic, Re(cr) = 0, one can distinguish two situations. In the excited case, 

a < 0, Tr pt^ converges to zero. In the dissipative case, a > 0, it converges to a finite value, which is not surprising 
since the system then reaches a thermal equilibrium [6]. For instance in the case of a bath of photons, one has 



Trpt' 



2n+l' 



(72) 



If the system is hyperbolic, the "Lyapunov exponent" w = Re(cr) ^ 0. The consequence is to shift the definition of 

2 

the above dichotomy. Indeed, Tr pt has a nonzero limit in the more restricted range a> u). 

The decoherence time, for the decay of Tr pt'^ , is in general inversely proportional to the coupling strength. However 
in the weak coupling limit \a\ <C w it is l/w in the hyperbolic case, whereas it is l/\a\ in the elliptic one. Hence 
the decoherence time defined by the linear entropy is here consistent with the positivity threshold. This is a strong 
support to the thesis that positive Lyapunov exponents accelerate decoherence [20] . 

The asymptotic formula (70) holds only for those cases where all the eigenvalues of M(— i) tend to infinity. Coun- 
terexamples are the elliptic case with a < 0, since the determinant then has a finite limit, and the hyperbolic case, 
for a < —CO. Moreover, though —oj < a < u; leads to a determinant which tends to infinity, one of the eigenvalues will 
indeed diverge whereas the other one, corresponding to the unstable direction, will have a finite limit. 
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VII. REVERSIBILITY OF THE SOLUTION 



Although decoherence is usually associated with a loss of information, which goes along with the growth of entropy, 
several techniques have been developed in quantum optics to recover the initial information after interaction of the 
system with the environment. In general these consist of reconstructing the quantum state of a lossy cavity by using 
mathematical inversion formulas [10] or directly by experimental processes [21]. We show here that the reversibility 
of the Wigner function results from the deconvolution of its evolution (24), or, even simpler, as a division (19) in the 
chord space. Indeed, one has 



Wt{0 = Woie-'^'R-tOGtiO, (73) 



which can easily be inverted as 



G,(e«'R,{) 

This is a generalization and a simplification of the inversion formula of [10] since the loss induced by beam splitter is 
a particular case of our general formalism. 



VIII. CONCLUSION 



The exact solution of the Markovian master equation for quadratic Hamiltonians and linear complex Lindblad 
operators has been derived in the form of a convolution for the Wigner function. This involves the classical evolution 
of the initial Wigner function for the closed system with a phase space Gaussian that is independent of this initial 
condition, while its width expands in time, depending only on the Hamiltonian and the Lindblad operators. This 
simple solution allows for the generalization of Diosi and Kiefer's proof that the Wigner function becomes positive 
within a definite time tp. Furthermore we support the much stronger statement that, unless the initial distribution 
is already a coherent state, the Wigner function must have negative regions before that time tp, which then does not 
depend on the initial state. In section V we have given the behaviour of tp for three basic types of motion of the 
quadratic Hamiltonian, namely the elliptic, the hyperbolic and the parabolic case, verifying that positivity is always 
reached exponentially fast. The positivity threshold is generally of the order of l/|a|, except in the weak coupling 
regime |a| <^ui of a. hyperbolic system, where it is of the order of 1/w. One should note that the threshold becomes 
independent on the Planck constant, if the Lindblad equation is appropriately scaled. 

The Fourier transform of the exact solution, the chord function, is even simpler. This is the product of two terms: 
one is just the nonhamiltonian classical evolution of the initial chord function with dissipation and the other is again a 
Gaussian, but with narrowing width. This leads to a simple formula for Tr fp as an average of a shrinking Gaussian, 
where the probability distribution used to calculate the mean is just the square modulus of the initial chord function. 
For long times, we use the normalization condition that the chord function is unity at the origin to derive simple rules 
for the asymptotic growth Tr for a general quadratic Hamiltonian. This decays exponentially fast in a time which 
is of the same order as the positivity threshold. This result is compatible with the arguments presented by Zurek 
and Paz [20] for exponential growth of linear entropy for chaotic systems. These arc classically characterized by local 
hyperbolicity, where the Lyapunov exponent describes the average effected by a typical orbit that approaches many 
hyperbolic points. In contrast, the hyperbolic quadratic Hamiltonian defines a linear classical motion, but both will 
exponentially stretch the Wigner function. Of course, in a chaotic system, the result must be analyzed more deeply 
since the phase space volume remains finite, which leads to a saturation of entropy even with no dissipation. 

We finally point out a very simple inversion formula which allows one to retrieve the initial state of the system, 
which seems a very transparent way to deal with the topics of quantum state reconstruction. 

We thank Luiz Davidovich, Ruinet Matos Filho and Fabricio Toscano for helpful discussions. We acknowledge 
financial support from Faperj, CNPq, Pronex and the Instituto do Milenio for Quantum Information. 

APPENDIX A: LINDBLAD EQUATION IN THE CHORD SPACE 

From the definition (4) of the chord transform, 

W{0 = exp (^-y A W(x), (Al) 
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applied to the derivatives of W and the product of W with q or p, we set the following transformation rules by using 
integration by parts: 

xW ^ __j 

dW ~ dW 
^._^.2W-(.^. (A2) 

By applying these rules on the following Poisson bracket, 

dW dW 
{if(x), W{x)} = 2(Hi2p + H22g)^(x) - 2(Hnp + Hi2g)^(x), (A3) 



we get 



that is 



dW dW 
{H{x),W{x)} ^ 2(Hi2Cp + H22C<,)^(x) - 2(HnCp + Hi2^J^(x), (A4) 



{H{x),W{x)}^[h{(,),W{0}- (A5) 

In the same way we get 

hence the final equation (10) for W, with the help of the last line of (A2). 

APPENDIX B: ZEROES OF THE HUSIMI FUNCTION 

Defining the complex variable ^(x) = {q + ip)/^/^, we may express the coherent state \z) [16] as 

|.)=e-NV2^A.^|n), (Bl) 



n=0 v/"' 

where |n) are the eigenstates of the harmonic oscillator. Hence, the coherent state representation of any pure state 
1'^) can be expressed as 

{z\i^) = e-l^l'/2F(0*), (B2) 

where 

is an entire function, known as the Bargmann function [22] . Since the Husimi function is the square modulus of the 
coherent state representation, we obtain 

Q{x) = \{z{x)m^=e-\^\'\F{z)\\ (B4) 

Thus the zeroes of the Husimi function coincide with the zeroes of an analytic function. Indeed, wc may define the 
state IV') by its Bargmann representation F{z). An important a priori restriction is that F{z) is at most of order 
r = 2. 
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To see this, recall that F{z) is of finite order if 

< el^l" (B5) 

for all sufBciently large \z\ and the order of this function is r = inf fj, for which (B5) holds. But if we use the fact that 
\{z\tlj)\^ < 1 in (B4), we obtain 

|F(^)| <el^l'/2 <el^l', (B6) 

so r < 2. 

We now make use of the following 

Theorem [23]: If F{z) is an entire function of finite order with no zeroes on the plane, then its order is necessarily 
an integer and F{z) = e^^'^''^ where Pr{z) is a polynomial of order r. 

In the case of the Bargmann function, then Pr{z) is at most of second order and hence given by (B4) must 

be a Gaussian in p and q if it represents a normalized function. 

Thus, only Gaussian Husimi functions have no zeroes in the phase plane. To a great extent, positions of the isolated 
zeroes of the Hiisimi function also restrict the class of admissible pure states through the factorization theorems of 
Weierstrass and Hadamard [23] . The characterization of the pure states as chaotic or regular by the pattern of zeroes 
has been extensively studied for the case where the phase space is a torus, because the restriction is then more severe 
[24]. 
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